Dark solitons in a superfluid Fermi gas 
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We investigate the behavior of dark solitons in a superfluid Fermi gas along the BCS-BEC crossover 
by solving the Bogoliubov - de Gennes equations and looking for real and odd solutions for the order 
parameter. We show that in the resonance unitary region, where the scattering length is large, the 
density profile of the soliton has a deep minimum, differently from what happens in the BCS regime. 
The superfluid gap is found to be significantly quenched by the presence of the soliton due to the 
occurrence of Andreev fermionic bound states localized near the nodal plane of the order parameter. 
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INTRODUCTION 

The interplay between effects of coherence and inter- 
action is one of the most interesting features exhibited 
by superfluids. It shows up, in particular, through the 
occurrence of topological excitations such as vortices and 
solitons. These have been the object of systematic inves- 
tigations in the case of ultracold Bose gases, where coher- 
ence is the result of Bose-Einstein condensation (BEC) 
and is associated with long range order in the one-body 
density matrix. The situation is even more interesting in 
the case of Fermi gases where coherence originates from 
the interaction between particles which, at low tempera- 
ture, brings the system into a superfluid phase character- 
ized by long range order in the two-body density matrix. 
For this reason, the consequences of coherence on mea- 
surable quantities, like the density distribution of the gas, 
are more indirect and subtle than for bosons. 

Fermi superfluid gases are now available experimen- 
tally in 40 K and 6 Li and various regimes along the BCS- 
BEC crossover can be explored by tuning the atomic 
s-wave scattering length a through the Feshbach reso- 
nances exhibited by these atoms. When a is small and 
negative the system is described by the BCS theory of 
superfluidity. Conversely when a is small and positive 
dimers of atoms of different spin are formed and the sys- 
tem behaves like a BEC of molecules. Near resonance 
the scattering length is much larger than the average in- 
tcrparticle distance and the system enters the so called 
unitary regime. Quantized vortices have been recently 
observed along the BCS-BEC crossover [l| and have been 
the object of several theoretical papers [3, 0, 0, H, E|- 

In this work we investigate another important class 
of nonlinear topological excitations: dark solitons. In 
three dimensions (3D) a dark soliton is characterized by 
a real order parameter which changes sign at a planar 
node (a point node in ID). In the BEC regime the soli- 
ton is a solution of the Gross-Pitaevskii (CP) equation 
for the orderparameter of the condensate with repulsive 
interaction [7(. In a uniform 3D system dark solitons 
are known to be unstable via the snake instability, i.e., 



a sinusoidal transverse oscillation of the planar node. In 
trapped gases, however, the instability timescale can be 
very long, so that solitons can indeed be observed ||. 

While in the BEC case the node of the order parame- 
ter causes a notch in the density distribution, in a BCS 
superfluid the density is almost unaffected by the pres- 
ence of the node [9( . The situation is similar to the one 
of vortices. As in that case, the behavior of the density 
along the BCS-BEC crossover is expected to be inter- 
esting and rather nontrivial, as a result of the delicate 
balance of coherence and nonlinear interactions. A ma- 
jor question concerns the behavior at unitary where no 
exact many-body theory is presently available. The prob- 
lem also shares useful analogies with the interference be- 
tween two expanding Fermi superfluid gases where the 
order parameter is expected to exhibit an oscillating be- 
havior with a change of sign, but no quantitative predic- 
tions are available concerning the behavior of the density. 



BOGOLIUBOV - DE GENNES THEORY 

We investigate the problem by using a mean-field the- 
ory for a 3D Fermi gas at zero temperature, based on the 
solution of the Bogoliubov - de Gennes (BdG) equations 
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where H = -h 2 X7 2 /2j\I + V cxt {r) - // is the single- 
particle grand-canonical Hamiltonian. For given chem- 
ical potential, fi, and order parameter of the super- 
fluid phase, A(r), these equations provide the spec- 
trum e v and quasiparticle amplitudes ^(r) and ^(r), 
which are required to satisfy the normalization rela- 
tion / d 3 r [u*,(r)u n (r) +v*,(t)v v (t)] = S v > n . The above 
equations must be solved together with the equations for 
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the order parameter and the density: 

A(r) = -flX^KM. (2) 
n(r) = 2^|^(r)| 2 . (3) 



This is done by means of an iterative procedure, which 
starts from a trial function A(r) and converges to the self- 
consistent solution of fl])-([3|). All sums in the equations 
are limited by an energy cutoff, < e n < E c . This cutoff 
is required in order to cure the ultraviolet divergences 
and is accompanied by a regularization of the interaction 
parameter g according to 111 ]: 
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Here a is the 3D s-wave scattering length characterizing 
the interaction between atoms of different spins, while 
Ef = h 2 kp/2M and Uf = (S^no) 1 ^ 3 are the Fermi en- 
ergy and momentum of a uniform ideal Fermi gas of den- 
sity no For sufficiently large values of E c the final 
results should not depend on the choice of the cutoff. 
It is worth stressing that the above equations reduce to 
the stationary GP equation for the order parameter of a 
condensate of molecules in the BEC limit , the inter- 
action between dimers being however given by the mean 
field value 2a instead of the exact value 0.6a 14{. Al- 



though approximate, this mean field theory is expected 
to give a comprehensive and reasonably accurate picture 
of the BCS-BEC crossover. 



RESULTS FOR A DARK SOLITON 

We look for solutions of CHJ)— dSJ) corresponding to a 
soliton at rest in the superfluid and in the absence of 
external potential (V ext = 0) The order parame- 

ter, A(z), is chosen to be a real and odd function of z, 
having a node in the xy-plane at z = 0. The density 
n(z) is an even function of z. The condition of reality 
implies the absence of currents, which are in general as- 
sociated with a z-dependence of the phase of the order 
parameter. The calculation is done in a finite box of size 
L x Lj_, where L is the size along z and all solutions 
are forced to vanish at the boundaries. If the box size is 
large enough the effects of the boundaries on the soliton 
are vanishingly small. Due to the translational symme- 
try in the transverse direction, the BdG eigcnfunctions 



can be written as u^(r) = u n ^ n± (z) e ± ' r± y/kp/L^ and 
v v (r) = v n>n± (z) e lk± ' r± \/%f/Lx, where r± = (x, y) and 
h± = (k x ,k y ). The transverse momentum is quantized 
according to the rule k x = 2irn x / L± and k y — 2Tm y /Lj_, 
with n± = (n Xl n y ) and n x , n y integers. We consider a 
gas in a box with L = 40k p 1 and Lj_ = 20k p 1 , and we 
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FIG. 1: (a) Order parameter and (b) density for a dark soli- 
ton with (kpa) -1 = and ±1. The value (fcfa)" 1 = (solid 
line) corresponds to unitarity, while (kpa)~ 1 = —1 (dashed) and 
(kpa) -1 = +1 (dot-dashed) are on the BCS and BEC side of the 
resonance, respectively. Both A(z) and n(z) are normalized to their 
asymptotic values far away from the soliton. 



also set E c = 50£f, which turns out to be large enough 
for a reasonable convergence. 

In Fig. Q] we show the order parameter (a) and the 
density (b) for a dark soliton with (fc^a)" 1 = —1, 
and +1. Both quantities are normalized to their asymp- 
totic values far away from the soliton, no = kp/3ir 2 and 
Ao. For Ao and n we find the values: Ao ~ 0.21ef and 
H ~ 0.96e F for (kpa)- 1 = -1 (BCS side); A ~ OJOep 
and \i ~ O.GIef for (kpa)^ 1 = (unitarity); Ao ~ 1.5ei? 
and m ^ -0.84e F for (^a)- 1 = +1 (BEC side). These 
values almost coincide with those of an infinite uniform 
system when calculated with the same value of the cut- 
off energy. This proves that the size of the box is large 
enough to neglect its effects on the calculations. They 
instead differ from the values calculated in a uniform gas 
in the limit of infinitely large E c . The difference is about 
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~ 10% on the BEC side and ~ 1% in the other cases. A 
larger discrepancy is of course obtained by using a lower 
cutoff energy, especially on the BEC side of the resonance 
where the formation of molecules with energy ~ h 2 /ma 2 
requires large values of E c to reach convergence. 

Figure [TJa shows the occurrence of a deep depletion of 
the density at unitarity (solid line) with a ~ 80% con- 
trast, comparable to the one of the BEC regime (dash- 
dotted). On the BCS side (dashed), conversely, the con- 
trast is only ~ 30% at (kpa)^ 1 = —1 and becomes expo- 
nentially small in the limit kp\a\ <C 1. These results are 
consistent with those obtained for the profile of a vortex 
core |,@|. 

In Fig. [Iji one notices that, as for the core of a vortex, 
the order parameter of the soliton in the BCS regime 
exhibits two length scales: a steep slope in a narrow re- 
gion of the order of kp l and a smoother slope in wider 
region of size £bcs = hvp/ Ao, wnere £bcs is the coher- 
ence length of the Cooper pairs and vf = \/2ep/M is the 



Fermi velocity. For (fcpa) 



-1 one has £bcs — Wk, 



We also find that the density exhibits oscillations with 
wavevector ~ 2kF- The same type of oscillations are 
found at the box boundaries and we checked that their 
shape does not depend on neither the box size nor the 
cutoff energy. One can thus safely identify them as 
Friedel oscillations. Approaching unitarity the coherence 
length £bcs decreases, eventually reducing to A;^ 1 . 

Let us focus now on the single-particle excitation spec- 
trum. In the uniform superfluid the spectrum is given 
by the well known result £buik(fc) = [(h 2 k 2 /2M — /i) 2 + 
A 2 ,] 1 / 2 This expression is plotted with solid lines 

in Fig. [2a-c) as a function of fcj_, for the three cases 
(fcpci) -1 = and ±1. The minimum of this func- 
tion defines the gap A gap for bulk excitations, which 
is A gap = A for fj, > and A gap = [/i 2 + A 2 ,] 1 / 2 for 
/i < 0. In the presence of the soliton, however, the eigen- 
values e n ,n_L °f t ne BdG equations exhibit a nontrivial 
feature: besides the continuum of bulk states with en- 
ergy above £buik one finds several states even below Ebulk- 
The energy of the lowest states for each k± is given in 
Fig. [2ja-c) (open circles). The corresponding quasi- 
particle amplitudes are localized near the soliton, as one 
can see in Fig.[2][d-f) where we plot the function \vo(z) \ 2 
of the state with k± = 0. The origin of these localized 
Andreev-like states [l?! resides in the fact that the en- 
ergy cost for creating a fermionic excitation near a node 
of the order parameter is reduced with respect to the 
bulk value. In the BCS limit the energy of the lowest 
bound state, also called minigap in the context of super- 
conductivity, is expected to be of the order of Aq/2ef- 
For a vortex in a BCS superfluid described by BdG equa- 
tions this result was proved in Ref. 18|, but the result 



is rather general and can be derived also when A(z) is 
a step function [13, |2l|- In the same limit, the eigenvec- 
tors Uo(z) an d vq(z) of the lowest bound state behave like 
~ cos(fc F z) exp(2:/^scs) and ~ sin(fc F z) exp(z/£ B cs), re- 
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FIG. 2: Left panels: Energy of the lowest single-particle states 
for fixed kj_ and for (a) (kpa)" 1 = —1, (b) and (c) +1. Circles 
correspond to the lowest bound states in presence of the dark soli- 
ton; solid lines correspond to the lowest energy states, £bulk> i n a 
uniform superfluid (see text). Right panels: |i>o(z)| 2 for the lowest 
bound state with kj_ = 0. The dashed line in (d) corresponds to 
the ansatz vq(z) ~ sin(fc^2;) exp(2:/^BCs) with Cbcs = lOfc^ 1 . 



spectively, as shown in Fig. [2ji. For the minigap at 
(k F a)~ l = -1 we find A m i mga p ~ 0.64Ag/2e_F. 

We now discuss the BEC limit, where the density 
profile should approach the analytic result obtained by 
solving the GP equation for a gas of bosons with mass 
rriB = 2M interacting with a scattering length cibb- 
The density of the bosonic dimers, rig = n/2, is given 
by 0: n B {z) = n B0 tanh 2 (z / \/2£bec ) , w ith n B o = 
n /2 = fc^/(67r 2 ), ^bec = fr/y/2m B gBBnB0 and g BB = 
4irh 2 a BB /m B . Analytic solutions can be found, in the 
same limit, also for the fermionic bound states. In fact, 
these states can be obtained by solving a Schrodinger 
equation for a fermionic impurity of mass mp = M in an 
inhomogeneous bosonic superfluid. By minimizing the 
total energy functional one derives the following equa- 
tion for the impurity wave function 'J: 



2m f 



dl + 9BFn B {z) 



(•5) 



where gsF = 2irh 2 a B F jm r is the dimer-atom coupling 
constants, while m r — m B mp / {m B + mp) — 2M/3 is 
the reduced mass. The exact values of the a BB and 
a B F scattering lengths are equal to 0.6a [lj| and 1.2a 
respectively, while the BCS mean field theory yields 
the Born approximation values 2a and 8/3a [23| . Equa- 
tion (|5|) corresponds to the Schrodinger equation for a 
particle in a modified Poschl- Teller potential. The so- 
lutions are analytic [24[ and include bound states. The 
energy of these bound states can be used to calculate the 
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case of Josephson currents through a potential barrier 
27] . In the opposite BEC limit, the binding of the un- 



FIG. 3: Dashed lines: density profile ub(z) /tibo and wave func- 
tion of the lowest fermionic bound state, ty()(z), in the BEC limit. 
Solid lines: density profile n(z)/no and lowest eigenvector 110(2) of 
the BdG equations for (fcj?a) _1 = +1. 



paired Fermi atoms on the solitonic plane shares an inter- 
esting analogy with the Andreev state of 3 He atoms on 
the free surface of superfluid 4 He |28[ . In the geometry of 
the soliton, the fermionic bound state can give rise to a 
2D Fermi gas embedded into a molecular Bose superfluid. 
This can be particularly interesting when the molecular 
BEC is made starting from a slightly imbalanced spin 
population. The residual gas of unpaired atoms can eas- 
ily fill the available states in the soliton, thus forming 
a polarized 2D Fermi gas. From the experimental view- 
point, it seems quite possible to observe solitons in the 
BCS-BEC crossover, for instance, by producing them in 
the molecular BEC phase with known techniques Q and 
then tuning the scattering length across the Feshbach 
resonance, similarly to what has already been done with 
vortices. A natural extension of this work is the study of 
solitons characterized by a complex order parameter and 
moving in the superfluid (grey solitons). 



minigap for fermionic excitations in the BEC regime. In 
fact the breaking of a pair causes the unbound motion of 
two fermions which will occupy the single particle states 
^0 with lowest energy. Using the mean field values for 
the scattering length, the energy of the lowest state turns 
out to be exactly one half of the bulk value [25| 
so that 

A gap - A minigap = 4irh 2 an /M . (6) 

Equation (|6|) provides a good quantitative estimate for 
the gap even for (fcpa) -1 = +1 where the above ex- 
pression gives the value A gap — A m ; n i gap = 0.26A gap to 
be compared with our numerical result ~ 0.27A gap . In 
Fig. [3] the analytic curve for the density ub (z) and the 
bound state wave function ^o(z) (dashed lines) are com- 
pared with n(z)/2 and uo(z) obtained from BdG equa- 
tions (solid lines). The agreement is rather good and the 
small differences are due to the fact that (fcpa)" 1 = +1 
is still relatively far from the true BEC limit and to the 
use of a finite cutoff energy E c . 

CONCLUSIONS 

In conclusion, we have theoretically investigated the 
behavior of dark solitons in a Fermi superfluid in the 
BCS-BEC crossover, showing that the soliton at unitarity 
has a large density contrast. We have also discussed the 
existence of bound states. Fermionic bound states have 
already been the object of investigations in the context of 
ID dark solitons in conducting polymers [5^|. They play 
an important role also in the physics of domain walls in 
the Fulde-Ferrell-Larkin-Ovchinnikov states (see, for in- 
stance, Ref. [ljl [20[ and references therein) and in the 
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